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ANALYSES OF THE SOLID EARTH AND OCEAN TIDAL 
PERTURBATIONS ON THE ORBITS OF THE 
GEOS-I AND GEOS-II SATELLITES 

T. L. Felsontreger 
J. G. Marsh 
R. W. Agreen 

ABSTRACT 

The luni-solar tidal perturbations in the inclination of the 
GEOS-I and GEOS-II satellite orbits have been analyzed 
for the solid Earth and ocean tide contributions. Precision 
reduced camera and TRANET Doppler observations spanning 
periods of over 600 days for each satellite were used to 
derive mean orbital elements. Perturbations duo to the 
earth's gravity field, solar radiation pressure, and atmo- 
spheric drag were modelled, and the resulting inclination 
residuals were analyzed for tidal effects. The amplitudes 
of the observed total tidal effects were about 1. 2 arc seconds 
(36 meters) in the inclination of GEOS-I and 4. 5 arc seconds 
(135 meters) for GEOS-II. The solid Earth tides were then 
modelled using the value for the 1^2 Love number of 0. 30 
available from Earth tide measurements, Earth rotation 
observations, and seismic data. The resulting inclination 


iil 



residuals were then analyzed for Kj, S 2 , and Pj ocean tide 
parameters. Since those ocean tidal constituents produce 
satellite orbital perturbations with distinctly different periods 
(ranging from 55 to 630 days), good separation of the con- 
stituent effects was possible. 

After the solution for these parameters, the inclination rms 
residuals were on the order of 0. 1 arc seconds (about 3 
meters). The satellite derived results for the Kj and 83 
ocean tide parameters are in good agreement with each other 
and with those deduced from existing surface data models; 
surface data models for the Pj tide were not available for 
comparison. The derived parameters consist of one second 
degree coefficient and an accompanying phase angle in a 
spherical harmonic expansion of the ocean tide potential 
for each tidal constituent — the results are as follows: 

Tidal Coefficient (cm. ) Phase 


Constituent 

GEOS-I 

GEOS-II 

GEOS-I 

GEOS-n 

S2 

Cj2= 1. 7±0. 5 

1 , 0+0.2 

e22=350°+21° 

62°+ll° 

Kl 

C2 i=8 . 8+0. 7 

5. 7+1. 6 

^ 21 = 15*^+ 4*^ 

334°±15° 

Pi 

021=5.0+1.1 

4.9+1. 1 

621 = 178°±15° 

127°+12° 
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ANALYSES OF THE SOLID EARTH AND OCEAN TIDAL 
PERTURBATIONS ON THE ORBITS OF THE 
GEOS-I AND GEOS-II SATELLITES 

INTRODUCTION 

The effects of the tidal deformations of the Eai’th on the orbits of close Earth 
satellites are readily observable as long period perturbations in the v-rbital 
elements of the satellites. This has been made possible by vastly improved 
tracking techniques and more accurate modelling of known perturbations such 
as the Earth's gravity field, lunar and solar gravity, solar radiation pressure, 
atmospheric drag, etc. Contributing to these tidal effects are the solid Earth, 
ocean, and (to a somewhat lesser extent) the atmospheric tides. However, 
since the same frequencios are present in the effects caused by these several 
contributions, any attempt to simultaneously recover parameters for tliem 
from satellite data is rather prohibitive. 

Previous work Involving analyses of tidal effects on satellite orbits has focussed 
on estimating values for the solid Earth tide Love number IC 2 and the associated 
phase lag. Early efforts wore undertaken by Kozai (1965, 1968) and Newton 
(1965, 1968). More precise work has been done recently by Anderle (1971), 
Douglas et al. (1974), and Smith et al. (1973). However, the results from 
these Investigations were satellite dependent and indicated a variance in k 2 
from a low of about 0. 22 to a higli of about 0. 31, which is contradictory to 
strong evidence pointing towards a value for 1^2 of about 0. 30. Attempts to 
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explain these discrepancies by the possible latitude dependence of the Love 
numbers (e.g. , Kaula (19C9)) have been disputed by Lambock (1974). The most 
feasible explanation for these anomalous results is the neglect of tiie ocean 
tidal effects, which Kaula (1962) first indicated could bo perceptible. Newton 
(1968) erroneously thought that these effects would almost cancel out when 
averaged over the whole Earth. More recently, Lambeck et al. (1974) have 
shown that neglect of the ocean tides could Introduce errors as large as 15 % 
in determinations of kg and sevet^'al degrees in the associated phase angle from 
analysis of satellite orbit data. Furthermore, as will be shown, in some 
instances even the atmospheric tidal perturbations cannot be neglected in 
determining solid Earth and ocean tidal parameters. 

In this paper, an attempt has been made to screen the solid Earth tide (and 
atmospheric tide) effects from satellite orbital data and analyze the result- 
ing data for ocean tidal parameters. The rationale for this approach is based 
on the assumption that at least the value for the second degree solid Earth 
tide Love number kg has been well established, while ocean tide models for 
the tidal constituents deemed to be Important are either in question or do not 
exist at all. Study has shown that long period tidal effects are most easily 
observable in the angle of inclination of close Earth satellites; therefore, the 
present work has been limited to analysis of this orbital element only. 


2 



TIDAL PERTURBATION EQUATIONS 


The ocean tide disturbing potential at a point (r, 0 , X ) exterior to the Earth, 
including the effects of ocean loading, can be expressed in a spherical harmonic 
expansion ns follows: 

io lo ■ 


• j^(a|^j 2 cos On ■'* On)cosgX 

+ (b|^g cos h d|^j 2 (1) 

where 

r = geocentric distance 

0 = latitude 

X = east longitude 

a,j = mean radius of the Earth 

g = acceleration of gravity 

P - mean density of sea water 

p^ - average density of tlie Earth 
K|^' = load deformation coefficient of degree k 

(?n = argument corresponding to a particular tidal constituent 

(Lambeck et al. , 1974; McClure, 1973) 

p 

•P|^ = Legendre functions 

^k£ * *^k£ ' ^k£ ' ‘^kfi ~ harmonic coefficients 
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Values for the load dofortnation coefficients are given by Longman (1066) and 
Farrell (1972). 


If we expand the potential of Equation 1 in terms of equatorial orbital elements 
of an artificial satellite, we obtain 




k+l k 

22 FkSp (I) G,jpq (e) 

p ==0 


Sin j^(k“2p) CO+ (k“2p+q) M-J* C (f2- Og) + + 

sin [(k“2p) w+ (k-2p+q) M+fi(f2-(?g)- + 

cos [(k-2p) (k-2p+q) M (H- 0 + 

+ Cj^jcos [(k-2p) + (k“2p+q) M •(• ^ ( S2- tfg) - c;^ - 


IcK ] I 

I k- 2 oven 

'«]) 

«m]| 

I k-£ odd 

‘^k«] I 


where 


( 2 ) 


a = semi-major axis of satellite's orbit 

I = angle of inclination of satellite's orbit 

e = eccentricity of satellite's orbit 

CO = argument of perigee of satellite's orbit 

M = mean anomaly of satellite's orbit 

n - longitude of ascending node of satellite's orbit 

Og = Greenwich mean sidereal time 
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“ inoikmtlon and occontrlclty functions defined by 
Knula (1966) 

®n “ ^ ” **n ^ corresponding to tidal constituent 

’f" * 

^IcK * ^kfi ° hariuonic coefficients 
®kie ' ^k8 ° angles 


Long period perturbations, i. e. , those having periods longer than one revolu- 
tion of the satellite, will occur only when the coefficient of M, k-2p+ q, is 0. 
In addition, perturbations having periods longer than a day will appear only 
when <ig disappears from the expression, I. o. , when £=■ n, which eliminates 
the terms multiplied by Cj^g . The principal perturbations occur when q = 0 
and k-2p = 0, so the resulting principal long period potential is 


f 1 


£=0 
k“£ oven 


(e) 


(«n- «5+ 


•f* *4* 

where Cj^g and Cj^g replace Cj^g and Cj^g , resp. For tlie semidiurnal 
tides, the most Important terms are those for which £ = 2 and k=2, 4, 6, ... . 
In the case of the diurnal tides, we have £ = 1 and k= 2,4,6 


Substituting the expression for in Equation 3 Into the Lagrangian equation 
of motion for I and integrating under the assumption that the argument is a 


5 


linonr function of time ylolds tho following perturbation equations for I, con- 
sidering only tho first harmonic coefficient for each tidal constituent: 


K) tide 
61 = - 

S2 tide 
51 = - 


P, tide 
61 = - 


9 


P 

■ (I+K2) 


cos I 


COS (S1+ 7T + C 2() 

10 


P<j) 

na^ (l-o^)2 

^Zl 

h 

9 


P 

(i+kJ) 

/ 3 

sin I 


Sin ( 2 J 2-2 x'+ 622) 

5 

'S' 

Pf!) 

\ " / 

na2{l-e^)2 

^22 

2 n - 2 \ ' 

9 


P 

, 1 

' % V 

cos I 


cos ( f2-2X’ + ff + c 

10 

S' 

P^ 



na2(l-o®)^ 

^2t 

n -2 x' 


( 4 ) 

(G) 

( 6 ) 


where 

n = mean motion of satellite 

^ ‘ecliptic longitude of tlie Sun 

Slmilor expressions can be derived for any of the other ocean tide constituents; 
however, as will be indicated Inter, strong evidence for observable effects 
from only tliese three constituents was found in tiie GEOS-I and GEOS-II data. 

In addition, since only the inclination data was considered, only one harmonic 
coefficient for each constituent could be recovered. This coeifioient would be 
biased by the effects of tho higher order coefficients (e. g. , C^j, . . . 
for the Kj tide), a fact which had to be taicen into account when comparing the 
results with existing surface data models. 
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computational procedures 


Derivation of ocean tidal coefficients from the analyses of GEOS-I and GEOS-II 
orbital data required fnvostlgntlon of tlio orbit evolution over several hundred 
days. The techniques and computer programs used to efficiently and accurately 
analyze observational data over orbital arcs of this length are discussed in the 
followlrig paragraphs. 

Orbit computations were carried out with the GEODYN Orbit and Geodetic 
Parameter Estimation Computer Program (Martin, 1972). The program em- 
ploys an 11th order Cowell numerical Integration procedure. Satellite disturb- 
ing forces modeled in the program Included the gravity field of the earth in 
the form of spherical harmonic coefficients, gravity anomalies, or surface 
densities, third-body gravitation, atmospheric drag, lunl-solar induced solld- 
Earth tides, and direct solar radiation pressure. For the present analyses 
the GSFC GEM-1 Earth Gravity Model (Lerch, ot al. , 1972) in the form of 
spherical harmonic coefficients was used. For the lunar and solar positions 
required in the computation of the third-body gravitational effects, the Jet 
Propulsion Laboratory (JPL) ephemeris tape was used. The atmospheric 
density was modeled by the Jacchia model atmosphere (Jacchia, 1965, 1968, 
1970, 1971). 

In order to investigate the long period orbital perturbations caused by the tides, 
a program designated ROAD (Wagner, et al. , 1974) was used. This program 
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effects of the solid Earth and ocean tides on the inclination. There is clearly 
evident an effect of about 1. 2 arc seconds (3G meters) in amplitude with a 
period of approximately 160 days, which is precisely the period of the Kj tidal 
effect, a luni-Bolar diurnal constituent. 

Figure 2 shows a similar result for the GEOS-n satellite. Here, a 651 day 
period (beginning March 1968) of both precision reduced camera and TBANET 
Doppler observations were processed In the form of two-day arcs to obtain 
mean orbital elements. The orbital elements derived from GEOS-II camera 
data were published b>> Douglas et al. (1973). The clearly visible effect has 
an amplitude of about 4. 5 arc seconds (135 meters). In addition to the solid 
Earth and ocean tidal perturbations, however, this effect also reflects a solar 
atmospheric tidal perturbation with an amplitude of 0. 1 arc seconds and a 
period of approximately 436 days. 

The following rationale for analyzing the data in Figures 1 and 2 was then 
adopted. Parameters for both the solid Earth and ocean tides cannot be re- 
covered simultaneously, since the same frequencies are Involved In each. 
However, the tidal effects seen in the data are essentially of second degree 
and the second degree solid Earth tide Love number has been fairly well 
established as a result of Earth tide measurements, Earth rotation observa- 
tions, and seismic data, while parameters for the ocean tide constituents are 
not too well known. Hence, the effects of the solid Earth tides were modelled 


preceding rage blank not filmed 
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using the value Ic2 ® .3 and a lag angle of 0 degrees, and the remaining 
Inclination residuals were analysed for ocean tidal effects. 

Figure 3 shows the GEOS-I inclination residuals after removal of the solid 
Earth tide effects. A frequency analysis performed on this data indicated the 
presence of perturbations having periods of 160, 85, and 56 days, which cor- 
respond to long period effects caused by the K|, P^, and S2 ocean tides, re- 
spectively. A subsequent least squares fit to this data involving the perturba- 
tion Equations 4, 5, and 6 resulted in the recovery of one harmonic coefficient 
and phase angle for each constituent (C21 and S2i for Kj and Pj, and C22 and €22 
for S2). Values for these coeificients and phase angles are given in Table 1. 

In addition, the amplitudes of the effects are given in meters (Figure 3) and 
arc seconds (Table 1). 

Figure 4 presents a similar set of Inclination residuals for GEOS-II after 
accounting for the solid Earth tidal effects. Also taken into account, however, 
is a substantial solar atmospheric tidal perturbation with a frequency about 
equal to that of the S2 ocean tide. The perturbation equation for this effect 
is very similar to Equation 5. Using the coefficient and phase angle for the S2 
atmospheric tide as reported in Chapman and Lindzen (1970) - C22 = . 352 mb, 
632 = 158° - it was found that this perturbation had an amplitude of slightly 
less than 0. 1 arc seconds. Again, a least squares fit to the data shown in 
Figure 4, for the Kj, Pj, and 83 ocean tides, yielded the results given in 
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Table 1. The periods and amplitudes (in meters) for these effects are shown 
on Figure 4 , while the amplitudes (in arc seconds) are presented in Table 1 . 

The rms of the residuals after the least squares fits shown In Figures 3 and 4 
were on the order of 0. 1 arc seconds for each satellite. These residuals rre 
shown in Figures 5 and 6. Frequency analyses indicate that there are no 
discernible significant periodicities left In these data. 

RESULTS AND CONCLUSIONS 

The coefficients and phase angles presented in Table 1 are in reasonable agree 
ment with each other, considering the differences in the amplitudes of the 
various effects from the two satelli„es and the fact that the effects of higher 
order harmonics have been "lumped" into the single coefficient. This latter 
point makes it somewhat incorrect to compare values for the same coefficient 
as derived from GEOS-I and GEOS-II. 

To illustrate this, let us write the perturbation equations for I considering the 
first two significant harmonic coefficients. They have the following form: 

K| tide 

61 = QCgi cos ( P + 7T + e2i) + IIC4J cos (P + w + e^j) ( 7 ) 

§2 tide 

61 = SCj2 sin (2 n -2 X '+ £22) + TC42 Sin(2 S 2 -2 \ ' + 642) (8) 

where 
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P 1 

f 

“ 10 

B - 


\ » , 

15 

p 

, / 

® ® \ 6 

16 ^ 


•( 1 +K 4 ) [ 

— ) 

9 

p 

, / 

tt ® Ns 

S= - g g 

P(S 

• ( 1 +K 2 ) 

— ) 

5 

p 

■H+k\) 

a® N 5 

T= g 

Pm 

-■) 

Equation 7 can be rewritten as 


® cos I 
na^ (1-e^)^ 


n 


na^ ( 1 -e^) * 
sin I 1 


n 


na® 2n-2 \' 


(9) 

( 10 ) 

( 11 ) 


61 = |^QC2iCos 621 RC41COS e^jj cos (12+ tt ) - |^QC2i sin 621 + RC41 sin e^jj 


sin ( J2 + 7T ) 


■( 13 ) 


or 

81 = VC21 cos ^21 ( ^"^ "■ ) “ V^21 ^21 ( i 2 + TT ) 

= VC21 cos (J 2 + TT + 62i)» ( 14 ) 

which is of the same form as Equation 4 . Similarly, for Equation 8, 

6 1= [sC22Cos £22+ TC42COS sin(2n-2X')+ [sC22 sin 622 + TC42 * 

• sin 642j cos (2J2-2x') 

= WC22COS 622 sin (2 J 2 -2 X ’) + WC22sin ^22 cos ( 2 f 2 - 2 X') ( 15 ) 

= WC22Sin (2 J 2-2 X ’+ ^22)* 
which is of the same form as Equation 5 . 
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Nowp 


V^2lC0S ^21 ® QC21COS Cji + RC^jCOS €41 ( 16 ) 

VC 2 J sin ^21 “ Q ^21 ^21 ^^41 ^4t 

W ^22 ^22 “ S ^22 ‘^22 ■*■ *^^42 ®°® ^42 ( 1 ®) 

W ^22 ®22 ~ ®^22 ^22 '^^42 ^42 ( 1 ®) 

Lambeck et al. (1974) have calculated values for the Cj^jg and from (1) 
the empirical solution of Dietrich (1944) for the Kj tide, and (2) from the 
numerical solution of Bogdanov and Magarik (1967) for the 83 tide. They are 
given in Table 2. 

If one computes values for Q, R, S, and T from Equations 9, 10, 11, and 12 
(using mean orbital elements), and uses the coefficients and phase angles in 
Table 2, the right hand sides of Equations 16, 17, 18, and 19 can be evaluated. 
Then, under the assumptions that V= Q and W= S, values for 621 , C 22 , ^ 2 i» 
and 622 flan be computed for each satellite. These are given in Table 3. 

These coefficients and phase angles provide a better basis for comparison 
with those in Table 1. (Unfortunately, no model for the Pj tide is available 
for comparison). There is reasonable agreement, except for the C 21 for 
GEOS-I and the C 22 for GEOS-II, Reasons for possible shortcomings in the 
assumptions made in this work are as follows: 
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(1) Inadequacy of the solid Earth tidal model. There Is some belief 
that the Love number IC 2 is not a constant for diurnal tides, but is 
frequency dependent due to core resonance and could Introduce 
errors on the order of 20% In the interpretation of tidal perturba- 
tions (McClure (1975)). 

(2) Attributing th j total effect, as seen in Figures 3 and 4, to ocean 
tides alone. 

(3) Inadequacy of the rjolar atmospheric tidal model for GEOS-H. 

(4) Neglect of liigher order coefficients in tlie solutions. 

(5) Inadequacy of the solar radiation pressure modelling. 

We believe, however, that these results indicate the feasibility of recovering 
global ocean tide parameters from satellite data. This capability will in- 
crease as data from ''nore satellites are studied, and as the effects are seen 
in other of the orbital elements (in particular, the longitude of the ascending 
node, for which the effects are greater than for the orbital inclination but are 
more difficult to isolate) . The main advantage in using satellite data for 
recovery of ocean tide parameters lies in the fact that the perturbation due to 
any tidal constituent appears as a long period effect with a frequency distinctly 
different from that of any other constituent, whether diurnal or semi-diurnal, 
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thereby making separation of the constituents quite easy. On the other hand, 
while surface tidal measurements can be used to separate diurnal from semi- 
diurnal constituents, it is difficult to differentiate between the various diurnal 
tor semi-diurnal) constituents unless observations over long periods of time 
are investigated. 
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Figure 4. Perturbation of the GEOS-n Orbit Caused by Ocean Tides 





Ocean Tide Parameters Derived from GEOS-I 
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Table 2 


Ocean Tide Parameters from Surface Measurements 


Source Constituent 

Dietrich (1944) Kj 

Bogdanov & MagarlU _ 

(1967) 2 


Coefficients (cm) 
C 21 " 3. 3| C 41 = 1. 2 

C 22 ~ 1" C^2 ~ 0» 2 


Phase Angles 
C 21 = 318°, ^41 =227° 

*=22 ~ 310 , £42 = 90 


Tidal 

Constituent 

Ki 

S2 


Table 3 

Combined Ocean Tide Parameters 


Coefftciep.ts (cm) Phase Angles 


GEOS-I 

GEOS-II 

GEOS-I 

GEOS-II 

^21 = 4. 0 

A 

C 24 ~ 5» 0 

C 21 = 351° 

%i = g' 

A 

A 

A 0 

A 

C 22 = 1* 7 

C 22 = 1.5 

£32 = 307 

^22 ~ 312 
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